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(57) Abstract: A method and system of improving sensitivity in the demodulation of a received signal over an arbitrary measurement 
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trajectory of the signal throughout the measurement time epoch (70); and the system comprising a receiver for receiving a direct 
sequence spread spectrum signal, the receiver comprising: an antenna ( 10) for receiving the direct sequence spread spectrum signal 
; a downconverter (40) for downconverting the received signal, producing a downconverted signal; an analog to digital converter(60) 
to convert the downconverted signal to a digital signal; a despreader (80) for despreading and coherently correlating the digital signal 
to a known signal, creating a despread signal; and a processor (70) for applying a Viterbi algorithm to the despread signal. 
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Field of the Invention 

The present invention relates to a method and system for increasing the 
sensitivity of the demodulation of direct-sequence spread spectrum modulated- 
signals and more particularly to the application of a Viterbi demodulation type 
algorithm to the output of a coherent correlation of a GPS C/A LI signal to improve 
the demodulation sensitivity. 

Background to the Invention 

Due to Federal Communications Commission (FCC) requirements for 
identifying mobile phone location, vendors are beginning to incorporate global 
positioning system (GPS) receivers into their handsets. This is an assisted form of 
GPS where the host base station provides key search parameter information to the 
handset such that the GPS satellite search and subsequent delay measurement 
becomes significantly less computationally extensive for the handset. In the code- 
division multiple access (CDMA) context, the messaging between the base station 
and the handset to support the assisted GPS (AGPS) is documented in "Position 
Determination Service Standards for Dual Mode Spread Spectrum Systems,* 
TIA/EIA/IS-80 1-1 /published by the Telecommunications Industry Association (TIA). 

In order for mobile location to be successfully based on AGPS, it is necessary 
that the GPS receiver in the handset unit have significantly higher sensitivity than is 
nominally required in a typical stand alone GPS unit. The reason is that the handset 
unit is usually oriented in positions or in locations that are not favorable for GPS 
signal reception. Hence a signal to noise ratio (SNR) sensitivity target is typically 
around 17 dB-Hz. Given that the GPS Coarse/Acquisition (C/A) L1 band which is 
typically used for mobile position location has a bandwidth of about 2 MHz, this 
implies that the input raw GPS, prior to processing, is more than 40 dB below the 
thermal noise floor. 

There is also a significant push to drive the cost of the overall handset down 
which limits the quantity of the processing which can be done on the received GPS 

« 

signals to extract them from the noise. This also limits the quality of the various 
radio frequency (RF) oscillator components used in the basket's receiver. 
Consequently, it can be expected that there will be some additional local oscillator 
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(LO) noise that needs to be compensated for. Notably, there is no closed automatic 
frequency correction (AFC) loop possible during the GPS measurement as the input 
SNR is too low. The lack of AFC, and the use of low cost, low power, RF LO 
components, implies that there will inevitably be a significant finite offset, drift, and 
5 some instability associated with the LO down conversion frequency. 

The C/A signal from a GPS SV (space vehicle) is a DS-SS (direct sequence 
spread spectrum) with a chip rate of 1 .023 Mcps. it is modulated as binary phase 
shift keying (BPSK) on a 1.574 GHz carrier. The GPS receiver correlates the 
received signal with a locally generated DS-SS code signal. In the AGPS scheme, 

10 traditional DS-SS correlation is also done. However, detailed information regarding 
the doppler shift of the SV GPS signal and code offset is available from the host base 
station (BS) which significantly reduces the search and detection effort. 
Nevertheless, the mobile GPS receiver is still required to determine the code delay 
and doppler to a finer resolution than that available from the host BS such that mobile 

15 location is possible. 

Based on standard assumptions regarding the noise in the GPS signal 
channel, the optimum receiver would correlate the signal in a coherent fashion over 
an integration time period that is sufficiently long to provide about 11 dB SNR at the 
correlator output. This will typically provide an adequate probability of detection with 

20 a reasonable false alarm rate. However, in the case of a mobile GPS receiver, due 
to the instability of the RF LO and the uncertainly of the SV doppler, the coherent 
integration epoch needs to be limited. Also valid GPS readings are still required 
even if the user does not hold the receiver steady. Hence, typically the coherent 
integration time is limited to 10 msec or less. As the available coherent integration 

25 epoch is not sufficient to obtain the sensitivity required, non-coherent summations of 
sequential coherent correlation integration outputs are used. However, non-coherent 
processing is a very inefficient means of further enhancing the SNR of a signal as it 
discards certain known statistical aspects of the signal. 

30 Summary of the invention 

The method and system of the present invention overcome the above by 
providing a way to demodulate the residual RF phase trajectory of a DS-SS 
dispreading and correlation, using quasi-coherent integration for a measurement 
epoch of arbitrary duration. The receiver in the present invention downconverts an 
35 incoming signal using a local oscillator. The receiver then despreads the signal. 
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In a perfect system the output from the despreader would be a DC signal. 
However, due to frequency and phase instabilities in the radio frequency local 
oscillator, and due to uncertainty of the GPS SV doppler, the phase trajectory is 
generally a time varying phasor. 
5 In order to provide increased sensitivity, a trellis is created with a fixed 

number of phase states evenly distributed between zero and 360°, The state 
trans/ilons in the trellis are limited based on known phase slew rate limitations. 
Probabilities for these state transitions are similarly based on these same factors. 

A Viterbi algorithm is then applied to the phase of the despread signal. 
10 Likelihood metrics are created for each path based on the phase of the despread 
signal as well as path probabilities. The system and method then uses the Viterbi 

■ 

algorithm to find the highest likelihood path for the phase trajectory through the trellis. 

The present invention therefore provides a method of improving sensitivity in 
the demodulation of a received signal over an . arbitrary measurement time epoch, 

15 said method comprising the steps of: correlating said received signal in a coherent 
fashion; and utilizing a Viterbi phase state keying trellis demodulation with a variable 
resolution of phase states over 0 to 360° to demodulate the radio frequency phase 
trajectory of said signal throughout the measurement time epoch. 

The present invention further provides a method of increasing sensitivity of a 

20 receiver for direct-sequence spread spectrum signals, said method comprising the 
steps of: downconverting the direct-sequence spread spectrum signal; despreading 
the downconverted signal creating a despread signal with a radio frequency phase 
trajectory; and applying a Viterbi algorithm to the despread signal to determine the 
radio frequency phase trajectory , whereby the algorithm has the steps of: defining a 

25 number of possible states, each state corresponding with a discrete phase and all the 
states together creating an even distribution of phases between zero and 360°; 
creating possible state transitions from each state; setting the probability of a given 
state transition occurring; defining a trellis with paths based on the possible state 
transitions at a number of discrete time intervals; taking a sample of the despread 

30 signal at a number of discrete time intervals; creating a likelihood metric for each 
path based on the phase trajectory of the despread signal at the time interval of the 
path; discarding all paths entering a state except for the path with the highest 
likelihood metric entering the state; and determining the most likely radio frequency 
phase trajectory based on the paths not discarded. 

35 The present invention still further provides a receiver for receiving a direct 

sequence spread spectrum signal, said receiver comprising: an antenna for receiving 
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the direct sequence spread spectrum signal; a downconverter for do wncon verting the 
received signal, producing a downconverted signal; an analog to digital converter to 
convert the downconverted signal to a digital signal; a despreader for despreading 
and coherently correlating the digital signal to a known signal, creating a despread 
5 signal; and a processor for applying a Viterbi algorithm to the despread signal. 

Brief Description of the Drawings 

Preferred embodiments of the present invention will now be described in 
greater detail and will be better understood when read in conjunction with the 
10 following drawings in which: 

Figure 1 is a schematic diagram of the receiver of the present invention; 

Figure 2 is an example trellis diagram for use with the method of the present 
invention; 

Figure 3 is an alternative trellis' diagram for use with the method of the 
15 present invention; 

Figure 4 is a graph illustrating an approximation of the probability distribution 
function of a phase sample as a periodic gaussian pulse on top of a constant 
function; and 

Figure 5 is a graph illustrating the losses in decibels verses the signal to noise 
20 ratio of the samples received in decibels. 

Detailed Description of the Drawings 

In this disclosure, a system and method are described that provide a practical 
way of performing an optimal correlation of the DS-SS GPS signal using the known 

25 statistical properties of the signal In particular, the method makes use of the 
generally known slew rate limitation of the LO instability and SV doppler. 

The main aspect of the method is a Maximum Likelihood Viterbi procedure 
similar to the processing used in trellis decoding of standard PSK communication 
signals. Although the method is applied exclusively to GPS signals in the following 

30 disclosure, one skilled in the art will realize that the method is also applicable for 
precision tracking of CDMA pilot signals for applications of high sensitivity mobile 
position location. 

Reference is now made to Figure 1 . The antenna 10 of a receiver 1 receives 
a direct sequence spread spectrum signal. The signal is then passed through band- 
35 pass filter 20 to isolate the signal. In a preferred embodiment, band-pass filter 20 is 
centred at 1 .574 Ghz, which is the frequency of the GPS C/A L1 carrier. 
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The output from band-pass filter 20 is then downconverted by multiplying it 
with the signal from the radio frequency local oscillator 30 in multiplier 40, producing 
a baseband signal. This baseband signal is then passed through low-pass filler 50 
and an analog to digital converter 60 to create a spread spectrum digitized signal. 
5 The spread spectrum digitized signal is then passed through a digital signal 

processor 70. Digital signal processor 70 is comprised of a despreader 80 and a 
processor 90. Despreader 80 is comprised of multiplier 84 and signal generator 83 
which generates coarse acquisition LI signal 82 which is multiplied with the spread 
spectrum digitized signal. 

10 In a perfect system, the output from despreader 80 would be a DC signal. 

However, due to frequency and phase instabilities in local oscillator 30, and due to 
uncertainty of the GPS SV doppier, the phase trajectory is generally a time varying 
phasor. The output of despreader 80 is thus processed by processor 90, which 
applies a Viterbi algorithm as described below. 

15 Consider the total GPS signal measurement time epoch, denoted by T^ as 

being divided into K equal length segments of T c each. It is assumed that the GPS 
correlator will coherently integrate for each of the segment lengths T c resulting in K 

correlation samples X=[x,,x 2 Xk] t . It is also assumed that T c is short enough such 

that uncertainties in GPS doppier, RF LO offset or frequency/phase instabilities of 

20 the RFLO are such that the resulting signal loss in the coherent correlation output 
samples in X is not significant. Typically T c will be somewhere between 1 to 10 
msec. 

X=[x 1f x 2 ,...,XK] is mapped into corresponding sample phases denoted by 
P=[Pi.P2,..mPk] t . These phases will be uniformly distributed between 0 and 360°. In 
25 order to implement a practical Viterbi algorithm, discrete phases are defined as being 
uniformly spaced such that: 

M = number of discrete phases 

s 0 ,s 1( ..-s w . 1 are the discrete phase states given by s n =(m/M) 360° 
Let k denote the k th state interval corresponding to the k m sample time or x k . 
30 Reference is now made to Figure 2, which shows an example trellis diagram for 
better illustration of the following description. 

Consider the m* state of the k* interval. If there are no phase dynamics due 
to the RF LO or doppier, then it will map into state m of the (k+1) m interval with a 
probability of 1. However, due to the uncertainty in the phase dynamics, the m th state 
35 can map into state (m-1) t m, (m+1) of the (k+1) m interval. Clearly, with these trellis 
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state transition probabilities, the modeled slew rate of the phase dynamics can be 
accurately limited. 

The trellis of Figure 2 is only one example. One skilled in the art will realize 
that other trellis state transition probabilities are possible, and that a probability trellis 
5 can be created to more accurately model these other probabilities. For example, in 
Figure 3, a system is illustrated in which the slew rate can be higher but a lower 
probability is assigned to higher slew rates. 

Intuitively such slew rate limiting lattice structures as described above 
effectively "low pass filter" the phase samples to reduce the influence of noise in 
10 correlation samples X. 

The number of phase states, H M\ is a function of the phase resolution 
required. If M is too small then there will be a small demodulation loss due to the 
discrete phase quantization. 

At each update, the standard Viterbi algorithm is used to eliminate the trellis 
15 paths of lower likelihood. (A good description of the Viterbi algorithm is found in 
"Digital Communications" by J. Proakis, McGraw Hill , 3 rd edition 1995.) Briefly, at 
the k* update a finite number of paths converge onto each of the M trellis phase 
states. For example, in the trellis structure of Figure 2, there are three candidate 
paths and in the trellis structure of Figure 3, there are five such paths. Each of these 
20 paths will have a likelihood metric associated with it. The path selected for each path 
is the one associated with the highest likelihood. Consequently, at each interval 
step, P paths converge to each trellis state and subsequently P new path segments 
are generated, as shown in the above examples for the (k+1) m trellis interval. 

The measured samples, p kl are the phases of the complex correlator outputs 
25 x k =Xik + ix^ as determined by: 

p k s atanfXqk/Xfc). 

This is an example only as there are a variety of numerical methods of 
extracting the phase from x k . In order to update the likelihood metric for each 
surviving path at the k* trellis interval, an estimate of probability of the received 

30 phase being p k conditioned on the phase state and the trellis path corresponding to 
the particular candidate path being considered is required. As the SNR of the 
correlator outputs, x k , is typically low (usually less than O dB), the probability 
distribution function (PDF) of the phase sample p k is not gaussian. Also the phase 
trajectory is dependent in a complicated statistical way on the candidate trellis path 

35 being considered. Consequently, finding the actual conditional probability is 
impractical. However, there are several assumptions and approximations that can be 
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reasonably made thai significantly simplify the likelihood probability metric. These 
are: 

1. Model the phase random process as first order Markov 

2. Approximate the PDF of p k as a periodic gaussian pulse on top of a 
5 constant function as illustrated in Figure 4. 

Given this approach, the conditional probability metric can be approximated 
as (pk-S m ) 2 where Sm represents the assigned phase of the m th state. The influence 
of the constant portion of P(Pk|S m ) is irrelevant as it is added to each candidate trellis 
path and therefore does not result in a path differentiating factor. Practical simulation 

10 experience has shown that this simple path likelihood metric is effective even though 
the approximation to the conditional probability of PtPntSm) is crude. 

Another very important factor is that the differentiating part of the likelihood 
metric when the above approximation is used becomes independent of the variance 
of the conditional probability of p*. This is a key factor in the robustness and 

15 practicality of the method. 

Example 

Below there is provided a simple MatLab simulation to demonstrate the 
effectiveness of the present method. 

20 In the simulation, an array of 2000 complex samples of x k were generated 

based on the assumption that the signal component was a constant amplitude 
exponential with a phase increment of 0.03 radians. Random noise was added to 
the samples of x k that was zero mean gaussian to both in-phase and quadrature 
phase components. The noise components were independent from sample to 

25 sample and also the noise in the quadrature and in-phase components were 
independent of each other. The variance of the noise components was selected to 
satisfy a given SNR of the correlation samples x n for the simulation. 

Tfie Viterbi demodulator was implemented as described above and tracked 
the phase excursion of the correlation samples x k . The accuracy of the demodulated 

30 phase was compared to the actual phase to determine a signal loss factor that was 
dependent on phase errors only. 

Figure 5 shows a plot of this loss factor as a function of the SNR of the 
correlation samples. As expected, the loss generally decreases as a function of 
increased SNR. At very low SNR levels, the phase tracking becomes so erratic that 

35 the Viterbi demodulation is not effective. 
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Figure 5 can be used to determine the number of samples required for a 
given SNR. Consider the samples at 0 dB SNR. As stated earlier, to satisfy the 
requirements of adequate detection and false alarm statistics, approximately 1 1 dB of 
SNR is required at the output of the overall correlation process. Hence if the 
5 samples of 0 dB SNR are generated after coherent correlation of T c , the number of 
samples required, K, is given by the formula: 

* 

10log(K) (Minimum Correlation Process Output SNR) - (Coherent 

Correlation SNR) + (coherence integration losses in dB) 

or, using Figure 5, 
10 10Iog(K) >- 1 1 - (0) + 1 = 12 dB. 

Soh/ing this, K=16. 

As indicated, the additional 1 dB added to the right hand side of the above 
formula is to offset the coherence integration losses represented in Figure 5. 
Comparing this with the traditional scheme which uses non-coherent summation of 
15 the correlator output samples, a minimum of K=26 is required. Hence, for a given 
GPS hypothesis dwell time, sensitivity can be significantly better with the Viterbi style 
demodulator as represented in this disclosure. 



20 Matlab simulation listing 

% viterbi demodulator for phase excursion" 



25 % number of state transitions 
K m 2000; 

% number of states 
M = 30; 

30 

% phase angle of each state 
state_angle - (2*D/M) * [0:M-1]\- 

% model the actual signal phase 
35 phi_sig = .03 * [1:K] ' ; 

% actual signal 
SNRdB = -6; 
SNR = 10" (0 . 1* SNRdB) ; 
40 u = exp{lj*phi_sig) ; 

x = randn(K,l) + 1 j*randn (K, 1) + sqrt (SNR) * u; 
phi_meas = angle (x) ; 

% metric for phase state 
45 % at each update there is only one surviving path for each 
state 
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Mt = zeros {M, 1 ) ; 

% paths, and initiate the k=l state; equal probability that we 
can start from 
5 % any state 

P = zeros (M,K) ; 
P<;,1) ^ [1:M] •; 
Pnew = P; 

10 % loop over the overall set of K state transitions 
• I = ones (M, 1 ) ; 

for k=2:K 

15 % determine the path metric of the new segment from state 

k~l to state k 

gl = phi_meas(k) * ones(M,l) - state_ angle; 

g2 = gl + 2*D*I; 

g3 = gl - 2 * Q*I; 
20 M_new = (min(abs( Egl,g2,g3] ■ ) ) ' ) .^2; 

% determine the new path for state m 
for m=l:M 

if m — 1 
25 iranl = M; 

mO = 1; 
mpl = 2; 
elseif m == M 
mml = M-l; 
30 mO = M; 

mpl = 1; 

else 

mml = m-l; 
mO = m; 
35 mpl = m+1; 

end 

[gl,i] s min( [Mt(mml) ,Mt(m0) ,Mt {mpl)] ) ; 
if i == 1 

40 Pnew(m,l:k) = (P(mml, 1 :k-l) ,m] ; 

Mt (m) = Mt(mml) + M_new(m) ; 
elseif i ==2 

Pnew(m,l;k) = [P(m0, 1 :k-l) ,m] ; 
Mt(m) = Mt(m0) + M_new(m); 

45 else 

Pnew(m,l:k) = [P(mpl, 1 :k~l) ,m] ; 
Mt{m) = Mt{mpl) + M_new{m) ; 

end 

end 

50 

P = Pnew; 



55 end % end of the K transition periods 
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figured) ;plot (P' ) ; 
hold on; 

[gl,i] = min(Mt) ; 
5 plot(P(i, :) , 'or' ) ; 
hold off; 

% unwrap P 
Puw = zeros (K, 1} ; 
10 Puw(l) = P(i, 1) ; 
offset = 0; 
for k=2:K 

Puw(k) = P(i,k) + offset; 
if Puw(k) - Puw(k-l) > 1 
15 offset = offset - M; 

Puw (k) = Puw(k) - M; 

end 

if Puw(k) - Puw(k-l) < -1 
offset = offset + M; 
20. Puw(k) b Puw(k) + M; 

end 

if abs(Puw(k) - Puw(k-D) > 1 

error (' discontinuity in Puw'); 

25 end 
end 



figure (2) ; 
30 plot (Puw); 

figure [ 3 ) ; 
plot(P(i, :> , % r« ) ; 
hold on; 

35 

plot (mod <phi_sig * M/ (2 * □),M), , g l ); 
hold off; 
grid on; 

40 

% calculate the integration loss 

Loss = 20 * loglO(abs(u f * exp (lj*state_angle (P( i, : ) ) ) ) /K) 



45 

The above-described embodiments of the present invention are meant to be 
illustrative of preferred embodiments and are not intended to limit the scope of the 
present invention. Also, various modifications, which would be readily apparent to 
one skilled in the art, are intended to be within the scope of the present invention. 
50 The only limitations to the scope of the present Invention are set forth in the following 

■ 

claims appended hereto. 
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I Claim: 

1 . A method of improving sensitivity in the demodulation of a received signal by 
a receiver over an arbitrary measurement time epoch, said method 
comprising the steps of: 

correlating said received signal in a, coherent fashion creating a 
correlation signal; and 

utilizing a Viterbi phase state keying trellis demodulation with a 
variable resolution of phase states over 0 to 360 Q to 
demodulate the radio frequency phase trajectory of said 
correlation signal throughout the measurement time epoch. 

2. The method of claim 1 , wherein said step of utilizing a Viterbi phase state 
keying trellis demodulation further comprises: 

breaking the measurement time epoch into time intervals; and 
creating a trellis of evenly distributed phase state nodes at each time 
interval, said creating step comprising: 

defining possible state transitions from each phase state node; 
creating state transition paths between phase state nodes in 
one time interval and phase state nodes in said one 
time interval; and 
assigning a transition probability to each state transition path. 

3. The method of claim 2, wherein the possible state transitions and the 
probability of the paths are assigned to reflect properties of said receiver. 

4. The method of claim 3 wherein the step of creating possible state transitions 
for each node is performed based on a known phase slew rate limitation of 
said receiver. 

5. The method of claim 4, wherein the known phase slew rate limitation is 
calculated from the instability of a radio frequency local oscillator in said 
receiver. 
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6. The method of claim 1, wherein the received signal is a global positioning 
system (GPS) coarse/acquisition LI signal generated by a space vehicle (SV). 

7. The method of claim 4, wherein the received signal is a global positioning 
system (GPS) coarse/acquisition U signal generated by a space vehicle (SV). 

8. The method of claim 1 wherein the received signal is a code-division multiple 
access (CDMA) pilot signal. 

9. The method of claim 7, wherein the GPS SV creates a doppler shift in the 
phase trajectory, and the known phase slew rate limitation is calculated from 
the uncertainty of the GPS SV doppler shift. 

10. The method of claim 7. wherein the known phase slew rate limitation is 
calculated from both the instability of a radio frequency local oscillator in said 
receiver and the uncertainty of the GPS SV doppler. 

11. The method of claim 2, wherein said step of utilizing a Viterbi phase state 
keying trellis demodulation further comprises: 

creating a likelihood metric for each trellis path based on a measured 
phase of the received signal and the probability for the path. 

* 

12. The method of claim 1 1 , wherein the likelihood metric is created based on an 
approximation of a probability distribution function of the phase of said 
correlation signal. 

13. The method of claim 12, wherein the approximation is to model the probability 
distribution function of the phase as a periodic gaussian pulse on top of a 
constant function. 

14. The method of claim 11, wherein the phase of the despread signal has a 
random process, and the likelihood metric is created based on an 
approximation of the random process of the phase of the despread signal. 

15. The method of claim 14, wherein the approximation is the modelling of the 
phase random process as a first order Markov. 

16. The method of claim 1 , wherein said receiver is a mobile receiver. 

12 



WO 2004/04.9663 



PCT/CA2002/001843 



17. A method of increasing the sensitivity of a receiver for direct-sequence spread 
spectrum signals, said method comprising the steps of 

downconverting the direct-sequence spread spectrum signal; 

despreading the downconverted signal creating a despread signal with a 

radio frequency phase trajectory; and 
applying a Viterbi algorithm to the despread signal to determine the radio 
frequency phase trajectory, whereby the algorithm: 
defines a number of possible states, each state corresponding with a 
discrete phase and all the states together creating an even 
distribution of phases between zero and 360°; 
creates possible state transitions for each state; 
sets the probability of a given state transition occurring; 
defines a trellis with paths based on the possible state transitions at a 

number of discrete time intervals; 
takes a sample of said despread signal at a number of discrete time 

intervals; 

creates a likelihood metric for each path based on the phase trajectory 
of said despread signal at the time interval of the path; 

discards all paths entering a state except for the path with the highest 
likelihood metric entering the state; and 

determines the most likely radio frequency phase trajectory based on 
the paths not discarded, 

18. The method of claim 17 wherein the direct-sequence spread spectrum signal 
is a global positioning system (GPS) coarse/acquisition (C/A) LI signal 
generated by a space vehicle (SV). 

19. The method of claim 17 wherein the direct-sequence spread spectrum signal 
is a code-division multiple access (CDMA) pilot signal 

20. The method of claim 17 wherein said receiver is a mobile receiver. 

21 . The method of claim 17 wherein the step of creating possible state transitions 
for each state is performed based on a known phase slew rate limitation of 
said receiver. 
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22. The method of claim 1 8 wherein the step of creating possible state transitions 
for each state is performed based on a known phase slew rate limitation of 
said receiver. 

23. The method of claim 21, wherein the known phase slew rate limitation is 
calculated from the instability of a radio frequency local oscillator in the 
receiver. 

24. The method of claim 22, wherein the GPS SV creates a doppler shift in the 
phase trajectory, and the known phase slew. rate limitation is calculated from 
the uncertainty of the GPS SV doppler shift. 

25. The method of claim 22, wherein the known phase stew rate limitation is 
calculated from both the instability of a radio frequency local oscillator in said 
receiver and the uncertainty of the GPS SV doppler. 

26. The method of claim 17, wherein the steps of setting the probability of a given 
state transition occurring and defining a trellis with paths based on the 
possible state transitions at each time interval are tailored to a particular 
receiver. 

27. The method of claim 18, wherein the likelihood metric is created based on an 
approximation of a probability distribution function of the phase of said 
despread signal. 

28. The method of claim 27, wherein the approximation is to model the probability 
distribution function of the phase as a periodic gaussian pulse on top of a 
constant function. 

29. The method of claim 17, wherein the phase of said despread signal has a 
random process, and the likelihood metric is created based on an 
approximation of the random process of the phase of said despread signal. 

30. The method of claim 29, wherein said approximation is the modelling of the 
phase random process as a first order Markov. 

31. A receiver for receiving a direct sequence spread spectrum signal, said 
receiver comprising: 



J4 



WO 2004/049663 



PCT/CA2002/001843 



an antenna for receiving the direct sequence spread spectrum signal ; 
a downconverter for downconverting the received signal, producing a 
downconverted signal; 

an analog to digital converter to convert the downconverted signal to a digital 

signal; 

a despreader for despreading and coherently correlating the digital signal to a 

known signal, creating a despread signal; and 

a processor for applying a Vlterbi algorithm to the despread signal. 

32. The receiver of claim 31, wherein the received signal is a global positioning 
system coarse/acquisition (C/A) LI signal generated by a space vehicle (SV). 

33. The receiver of claim 31 , wherein the received signal is a CDMA pilot signal. 

34. The receiver of claim 31 , wherein the receiver is a mobile receiver. 

35. The receiver of claim 31 , wherein the known signal is a GPS C/A L1 signal. 

# 

36. The receiver of claim 31 , wherein the processor 

breaks the despread signal into time intervals; arid 

creates a trellis of evenly distributed phase state nodes at each time 

interval by: 

defining possible state transitions from each phase state node; 
creating paths between phase states node in one time interval 

and phase state nodes in the time interval; and 
assigning a probability to each path. 

37. The receiver of claim 36, wherein the possible state transitions and the 
probability of the paths are assigned to reflect properties of said receiver. 

38. The receiver of claim 37, wherein the possible state transitions for each node 
are based on a known phase slew rate limitation of said receiver, 

39. The receiver of claim 38, wherein the known phase slew rate limitation is 
calculated from the instability of a radio frequency local oscillator in said 
receiver. 
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40. The receiver of claim 32, wherein the GPS SV creates a doppler shift in a 
phase trajectory of said received signal, and the known phase slew rate 
limitation is calculated from the uncertainty of said GPS SV doppler shift. 

41 . The receiver of claim 32, wherein the possible state transitions for each node 
are based on a known phase slew rate limitation of said receiver. 

42. The receiver of claim 41, wherein the known phase slew rate limitation is 
calculated from both the instability of a radio frequency local oscillator in the 
receiver and the uncertainty of the GPS SV doppler. 

43 The receiver of claim 36, wherein the processor further creates a likelihood 
metre for each path based on a measured phase of said received signal and 
the probability for the path. 

44. The method of claim 43, wherein the likelihood metre is created based on an 
approximation of a probability distribution function of the phase of said 
despread signal. 

45. The method of claim 44, wherein said approximation is to model the 
probability distribution function of the phase as a periodic gaussian pulse on 
top of a constant function. 

46. The method of claim 36, wherein the phase of the despread signal has a 
random process, and the likelihood metre is created based on an 
approximation of the random process of the phase of the despread signal. 

47. The method of claim 46, wherein said approximation is the modelling of the 
phase random process as a first order Markov. 
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